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ABSTRACT 

The baryon content of high-density regions in the universe is relevant to two critical unanswered 
questions: the workings of nurture effects on galaxies and the whereabouts of the missing baryons. In 
this paper, we analyze the distribution of dark matter and semianalytical galaxies in the Millennium 
Simulation to investigate these problems. Applying the same density field reconstruction schemes as 
used for the overall matter distribution to the matter locked in halos we study the mass contribution 
of halos to the total mass budget at various background field densities, i.e., the conditional halo mass 
function. In this context, we present a simple fitting formula for the cumulative mass function accurate 
to < 5% for halo masses between 10 10 and 10 15 /i _1 M Q . We find that in dense environments the halo 
mass function becomes top heavy and present corresponding fitting formulae for different redshifts. We 
demonstrate that the major fraction of matter in high-density fields is associated with galaxy groups. 
Since current X-ray surveys are able to nearly recover the universal baryon fraction within groups, 
our results indicate that the major part of the so-far undetected warm-hot intergalactic medium 
resides in low-density regions. Similarly, we show that the differences in galaxy mass functions with 
environment seen in observed and simulated data stem predominantly from differences in the mass 
distribution of halos. In particular, the hump in the galaxy mass function is associated with the central 
group galaxies, and the bimodality observed in the galaxy mass function is therefore interpreted as 
that of central galaxies versus satellites. 



Subject headings: cosmology: large-scale structure of universe 
numerical 



galaxies: groups: general — method: 



1. INTRODUCTION 

The environmental dependence of galaxy properties, 
such as broadband color, star formation, an d stellar 
mass , is a well-known effect in the local universe (jDresslerl 
1980). In this context environment means an estimate 
of the smoothed density filed at a given location. Of 
particular interest for the present study is the depen- 
dence of the galaxy stellar mass function (GMF) on envi- 
ronment. Recent comprehensive galaxy redshift surveys 
have led t o inte nsive studies in this field. For instance, 
I Mo et akl (|2004f ) model the dependence of the luminosity 
function on the large-scale environment based on mock 
catalogs for the t he two-degree Field Gal axy Redshift 
Survey f2dFGRS: IColless et al.ll200l and iBaldrv elaTl 
(|2006f ) investigate the GMF as a function of environ- 
ment in the nearby uni verse based on t he Sloan Digi- 
tal Sky Survey (SDSS; lYork et all l200"oh . To uncover 
evolutionary effects, surveys spanni ng a larger redshift 
range have been used: the study bv lBundv et al.l (|2006T ) 
is based o n the DEEP2 Galaxy Redshift Survey (0.4 < 
z < 1.4): iPannella et all (12001 use the COSMOS sur- 
vey; iBolzonella et alj (|2009t) "employ the zCOSMOS sur- 
vey in the redshift range < z < 1; and lScodeggio et ah! 
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(|2009[ ) investigate the environment dependence of the 
GMF based on the VVDS survey covering a redshift 
range of 0.2 < z < 1.4. 

On intermediate scales, ~ /i _1 Mpc, groups and clus- 
ters of galaxies themselves provide a definition of en- 
vironment. At low redshift, iBalogh et al.l (|2001f ) using 
Two Micron All Sky Survey (2MASS) and Las Campanas 
Redshift Survey (LCRS) data, separate different environ- 
ments as field, groups, and clusters, finding that galaxy 
luminosity and mass functions depend on both, galaxy 
type (with steeper functions for emission line galaxies) 
and mass of the group (with more massive and brighter 
objects more common in clusters), mainly as a conse- 
quence of the different cont ribution of passive galaxies 
(see also lHansen et al.ll2009|) . 

The environmental dependence of properties of gas 
is linked to the studies of warm-hot intergalac- 
tic media (WHIM). Hydrodynamical simulations by 
ICen fc Ostrikerl (|1999f) show that the average temper- 
ature of baryons is an increasing function of time, with 
most of the baryons at the present time having a tem- 
perature in the range of 10 5 - 10 7 K. The detection 
of this warm-hot gas poses an observational challenge. 
While according to their census more than one-half of 
the normal matter is yet to be detected, it is not clear 
which methods have to be used. In a later study 
ICen fc Ostrikerl (|2006l ) report that the gas density of the 
warm-hot intergalactic medium is broadly peaked at a 
density three to seven times the critical density, however 
their dark matter mass resolution was only moderate and 
a question of assigning baryons to groups has not been 
addressed. 

All the above studies would greatly benefit from a 
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quantification of the contribution of galaxy groups to 
the density fields as a function of background density. 
The signific ance of this question h as been recognized 
before fe.g.. ISheth fe Tormenlll999l ). but a detailed an- 
swer has never been provided. Simulation work has in- 
stead concentrated on questions regarding halo assembly 
bias and, in particular, on changes in the properties of 
galaxies within the halos of similar mass but residing 
in different environments (ILemson fe Kauffman"n] H99l 
iGao et al.ll2005t iCroton et al.ll2007ft . " 

This paper is organized as follows. In Section [2] we re- 
view the Millennium Simulation, the semi-analytic mod- 
eling of galaxies and the determination of the background 
density. Section [3] examines the dependence of the dark 
matter halo mass fractions on environment. In Section [4] 
halo mass functions in different environments are dis- 
cussed. Section[S]focuses on the dependence of the GMFs 
on environment. We present a short summary in Sec- 
tion [B] 

2. DATA 

This analysis is based on the pu blicly available 
Millennium simulation run (MS; ISpringell [20051: 
ILemson fe Springe! |2006|) . In the first part of this sec- 
tion, we review the MS and the scmianalytic galaxy mod- 
eling. The second part describes the determination of 
background densities and halo mass fractions. 

2.1. Millennium simulation 

The MS adopts concordance values for the parame- 
ters of a flat A cold dark matter cosmological model, 
f^dm = 0.205 and Ob = 0.045 for the current densities in 
cold dark matter and baryons, h = 0.73 for the present 
dimensionless value of the Hubble constant, as = 0.9 
for the rms linear mass fluctuation in a sphere of radius 
8 h~ 1 Mpc extrapolated to z = 0, and n = 1 for the slope 
of the primordial fluctuation spectrum. The simulation 
follows 2160 3 dark matter particles from z = 127 to the 
present day within a cubic region 500 ft,~ 1 Mpc on a side. 
The resulting individual particle mass is 8.6x 10 8 ft _1 M0. 
The gravitational force has a Plummer-cquivalent co- 
moving softening of 5 /t~ 1 kpc. The Tree-PM A^-body 
code GADGET2 (jSpringel et al.ll2005[) has been used to 
carry out the simulation and the full data are stored 64 
times spaced approximately equally in the logarithm of 
the expansion factor. 

The halos are found by a two-step procedure. In 
the first step, all collapsed halos with at least 20 par- 
ticles are identified using a standard friends-of-friends 
(FoF) group-finder with linking parameter b — 0.2. 
These objects will be referred to as FoF-halos. Then, 
post-pr ocessing with the su bstructure algorithm SUB- 
FIND (jSpringel et al.l 1200 ll ) subdivides each FoF halo 
into a set of self-bound sub-halos. Based on their 
assembly histories, individual sub-halos are populated 
with galaxies by a semi-analytic prescription and 
various observable quantities are generated. For a 
detailed description of the semi-analytic galaxy cat- 
alog we refer the r eader to ICroton et al.l (|2006[) and 
iDe Lucia fe Blaizotl (|2007l) . The stellar mass functions in 
this study are based on the DeLucia2006a_S DSS2 MASS 
catalog (|Lemson fc Virgo Consortium! 120061 

http : // www . g-vo . org/Millennium) . 

The following analysis is based on FoF halos which 



hereafter are addressed simply as halos. In principle, 
one also could use different halo definitions like those de- 
rived from spherical top-hat overdensity criteria or gravi- 
tational sclf-boundcdness. The former would case a com- 
parison with observations and the latter would not suffer 
from the bridging problem inherent to the FoF approach. 
However, according to various tests which we have per- 
formed, the results based on the various halo identifi- 
cation schemes show little difference. The advantage of 
using FoF halos lies in the simplicity of the approach 
which makes it a very common tool for the analysis of 
A^-body simulations. 

2.2. Background Density and Halo/Galaxy Mass 
Fraction 

The MS database also provides information on the 
global density field besides describing individual halo 
properties. Here, we use the densities, pcic, which are 
based on a Counts in Cell (CiC) approach using cubic 
cells of ~ 2 /i _1 Mpc on a side (jHocknev fc Eastwood! 
fl98l and PGio (hereafter background density), which 
are derived from the former by smoothing them with 
a 10 /i _1 Mpc Gaussian kernel. In this study, we utilize 
Pcic to compute the total amount of matter in a given 
volume and calibrate the use of the 10/i _1 Mpc back- 
ground density fields, pcio- We assume that pcio repre- 
sents the linear density field. 

To determine the total mass of halos per cell, we re- 
implement the CiC approach only taking into account 
the particles belonging to FoF halos above a given mass 
limit. If a halo crosses the boundary of a cell, only the 
mass of the halo inside the cell is counted. This quantity 
is used to compute the halo mass fractions in Scction [3.2l 
For the determination of the halo and the GMFs (Sec- 
tion[l]and Section[5]), we attribute the total mass of halos 
or galaxies to the cell within which their centers are lo- 
cated. Eventually, each cell is assigned a total mass in 
halos, a total mass in stars (galaxies), a total mass, and 
a smoothed background density (which after multiplica- 
tion with the cell volume corresponds to a mass as well). 

The upper left and the two lower panels of Figure Q] 
display a slice of the density contrast, (5qio = (pgio — 
(p))/(p) scaled by the variance of the matter fluctua- 
tions within a volume corresponding to the Gaussian fil- 
ter, (TgiOi at redshifts z = 0, z = 1, and z = 3. Typ- 
ical patterns of the large-scale density field are appar- 
ent, such as roughly spherical high density regions, fila- 
ments, and voids. According to linear theory, 5 /a does 
not change with time. Indeed, the scaled density contrast 
is very similar in shape and amplitude at the redshifts 
we consider. Figure [2] shows the differential mass distri- 
butions as functions of the scaled density contrast. Red 
curves show the lognormal fits to the di stribution (cf., 
iColes fc JonellllHlNevrinck et al.|[2009h . The approx- 
imate invariability of the scaled contrast reaffirms the 
presumption of linearity of <5gio- 

The top right panel of Figure [T] shows the halo distribu- 
tion within a slice of 10 ft, _1 Mpc thickness at z = at the 
same location as slice used for the contrasts. The under- 
laid gray scale image is a replication of the contrast on 
the left. As expected, the halo distribution follows the 
pattern of the background density field. In the following, 
we measure the dependence of the fraction of mass cap- 
tured in halos or galaxies as a function of the background 
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Fig. 1. — Upper left panel: dark matter density field at z = 0, smoothed using a Gaussian filter with a smoothing scale of 10 h _1 Mpc, 
5qioi an d scaled by the average rms fluctuations within the corresponding volume, ctqiq. Upper right panel: density field as shown on the 
left in gray scale. The points represent halos above 10 12 /i —1 Mq within a slice of 10h -1 Mpc thickness. Colors (from green to red) and 
sizes (from small to large) correspond to the logarithmic mass of the halos. As expected, the halo distribution closely follows the density 
field. Lower panels: smoothed and scaled density fields, as the upper left panel, but for the redshifts 2 = 1 and z = 3. According to linear 
theory 8/a does not change with time. 



density. 



3. HALO MASS FRACTIONS 



The first part of this section notes some general fea- 
tures associated with the characterization of the density 
field. In the subsequent paragraph, we examine the de- 
pendence of the halo mass fraction on environment. The 
halo mass fraction is defined as the mass locked in ha- 
los divided by the total mass in a given volume. As 
stated before, the environment is quantified based on the 
smoothed density within this volume. The halo mass 
fraction is closely related to the cumulative halo mass 
function which will be discussed in the subsequent sec- 
tion. 



3.1. Mass within Isodensity Surfaces 

The upper panel in Figure [3] shows the ratio between 
the average CiC-density of all cells located in regions of a 
given density contrast, (5gio, as a function of the density 
contrast itself. Line styles correspond to different red- 
shifts. At a contrast of 5gio ~ 0.2, the mass within a cell 
approximately corresponds to the value of pcio multi- 
plied by the cell volume. For higher contrasts, the actual 
mass deposited in the cell is larger than that deduced 
from the smoothed density field. The opposite is true for 
contrasts below 0.2. This behavior simply reflects that 
extremes are leveled by the smoothing procedure. 

More important for the subsequent analysis is the 
lower panel in Figure [3l It displays the mean value of 
Pcic/pgio within a volume confined by the isodensity 
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Fig. 2. — Differential mass distribution as a function of scaled 
density contrast. The density is computed based on a 10 h~ x Mpc 
Gaussian smoothing kernel, and the scaling factor, <jq w , is the 
mass variance within a corresponding volume. Solid, dotted, and 
dashed lines indicate the redshifts z = 0, 1 and 3, respectively. The 
red curves show the lognormal fits to the distribution. 
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Fig. 3. — Upper panel: mean CiC density as a function of 
smoothed density. Lower panel: mean density within isodensity 
contours divided by that isodensity. 

surface, pciOi as a function of pew or it's equivalent 
5qio- The ratio, pcic / Pcio, gives the average factor with 
which the confining density has to be multiplied to get 
the true mean density or the true mass within the en- 
closed volume. For different redshifts, indicated by dif- 
ferent line styles, one finds slightly different ratios, but 
the overall behavior is similar. For density contrasts be- 
tween and 1.5, which are of interest here, one obtains 
1 ^ (°cic//°gio is 3. This factor has to be accounted for 
when the total mass within a given volume is inferred 
from the confining isodensity surface. Loosely speaking, 



a factor of 2 has to be multiplied to the value of the con- 
fining surface density, PgiOj to recover the mass inside. 

3.2. The dependence of the halo mass fraction on 
environment 

The three panels in Figure [4] show the halo mass frac- 
tions for the redshifts z = 0, z = 1, and z = 3. The halo 
mass fraction is defined as the fraction of mass locked in 
halos (above a given limiting mass) and the total mass 
within a volume confined by isodensity surfaces with a 
density contrast <5gio as a function of that contrast. Line 
styles refer to halo mass limits as indicated. For each red- 
shift, these mass limits correspond to fixed values of the 
equivalent peak height, v(M,z) = 5 c (z)/a(M, z), where 
cr(M, z) is the rms linear overdensity within a sphere con- 
taining the mass M at rcdshift z, and S c (z) is the linear 
overdensity threshold for collapse at that rcdshift. The 
usage of v instead of the actual halo mass should remove 
much of the cosmology dependence of our results and 
will be most useful for the discussion of the halo mass 
functions below. The mass limits for the left and mid- 
dle panels correspond to v « 0.9, 1.0, 1.2, 1.4, 1.7, and 
2.1. In the z = 3 panel (on the right) mass resolution 
only allows us to show the graphs for the three highest v 
values, 1.4, 1.7, and 2.1. 

Figure [5] shows the fraction of mass in halos with re- 
spect to the total mass for confining density contrasts of 
(5gio = 0.5, 1.0 and 1.5 as a function of the mass limit 
of the halos. It is formally identical to the conditional 
cumulative mass function, F{v), which is discussed in 
the next section. However, here it is computed based on 
the particle distributions of halos. Therefore, the mass 
of halos which transgress cell boundaries is accurately 
split between the cells. The analysis in the next sec- 
tion is somewhat more coarse in the sense that it at- 
tributes the total mass of a halo to only one cell, namely 
the cell where the center of the halo is located. This 
can introduce some bias, similar in origin to the devia- 
tion of surface and mean enclosed density discussed in 
13.11 However, the comparison between the two methods 
yields acceptable agreement which implies confidence in 
the computation of the mass functions presented below. 

3.3. Where are the Warm-Hot Baryons? 

Our findings have consequences for the detection of 
warm- hot intergalactic baryons (WHIM). Cosmological 
simulations predict that some 50% of all the baryons lo- 
cally appear in the f orm of gas with temp eratures be- 
tween 10 5 and 10 7 K (jCen fc Ostrikerlll999t) . The main 
process responsible for heating the baryons to such tem- 
peratures is large scale structure formation and thus the 
missing component is thought to be associated with high 
densities. Observational detection of this component be- 
came subject of a number of studies, and yet a full ac- 
count of warm-hot bary ons has not been r eached. 

In a subsequent study iDave et"ail (|2001[ ) have argued 
that the major fraction of WHIM is located outside 
galaxy groups. Using hydrodynamical simulations they 
have shown that the total fraction of WHIM associated 
with galaxy groups is between 10% and 25%, depend- 
ing on the group mass cut and that the majority of the 
WHIM resides at mean densities of 10 < 6 < 200. In this 
work, no attempt has been undertaken to study the con- 
tribution of baryons in galaxy groups to the WHIM as a 
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Fig. 4. — Halo to total mass fraction as a function of background density contrast , 5gio- As indicated, the line styles correspond to 
different lower mass limits for the halos used to compute to total halo mass budget. For the left and middle panels, these masses correspond 
to the equivalent peak heights v = 0.9, 1.0, 1.2, 1.4, 1.7, and 2.1. In the 2 = 3 panel on the right only the limiting masses corresponding 
to v = 1.4, 1.7, and 2.1 are shown. At that redshift, i/-valucs below 1.4 correspond to halos with less than ~ 20 particles which are not 
resolved. 

side groups indica te that the universal baryon fraction is 
nearly recovered (iKravtsov et "all 120051 ISun et all [20091: 
iGiodini et all l2009t iMcGaugh et al.ll2010D . and poten- 
tially can be fully resolved by accounting the baryons 
near the virial radius. Thus, on the supposition that gas 
follows the dark matter our findings indicate that 20%- 
30% of the WHIM is located in groups if no background 
density constraint is imposed, whereas 60%-70% of the 
WHIM is residing in galaxy groups at high background 
de nsities. The forme r agrees well with the figures quoted 
in lDave et all (|2001| ). 

The only observational result p ublished on the global 
fractio n of mass in groups, by iReiprich &: Bohringerl 
(|2002t) . uses ROS AT All Sky Survey data, and quoted 
^duster = 0.012 (or 5% contribution of clusters to the 
total mass budget) for masses exceeding 10 14 Mq/i _1 . 
This measurement compares well with our value for 
M > 2 x 10 14 /i _1 M Q , which corresponds to the same v 
when a difference in the a% value between observed uni- 
verse and the Millennium run is taken into account. This 
account will likely improve soon since current X-ray sur- 
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Fig. 5. — Conditional cumulative halos mass function, F(y), 
i.e., fraction of mass locked in halos with respect to the total mass 
within a region confined by the contrasts, <5gio = 0.5, 1.0, and 1.5 
as a function of halo mass given parameterized by the equivalent 
peak height v = 8 c (z)/cr(M, z). The mass scale at the top indicates 
the equivalent masses for 2 = 0. Solid, dotted, and dashed lines 
correspond to redshifts z = 0, 1, and 3, respectively. At z = 3 <5gio 
does not reach 1. 

function of environment. Here we aim to determine the 
fraction of WHIM located in groups at high background 
densities. 

Figure U suggests a dominant contribution of groups to 
the matter budget at high overdensities, resolving 60% 
of total mass in groups with mass exceeding 2 x 10 13 M Q 
and 70% of total mass in groups with mass exceeding 
4 x 10 12 M Q . Without any background density restriction 
these groups account for 20%-30% of matter. 

It is important to note that no gas component has been 
included in the Millennium Simulation. To derive con- 
clusions on the WHIM, we rely on a model describing 
the distribution of gas relative to the overall matter dis- 
tribution which is equivalent to dark matter distribu- 
tion in the current context. Here, we assume for sim- 
plicity that gas follows dark matter. Hydrodynamical 
simulations show that this assumption is well justified 
for local densities below < 10 4 ti mes the cosmic mean 
density (|Faltenbacher et all 120071 cf.,[). Furthermore, 
observational and theoretical accounts for baryons in- 



veys can ac cess masses below 10 13 M Q (|Finoguenov et all 



120071 12009D . Also, detections of WHIM in X-ray emission 
regions agree well with o ur findings at high b ackground 
densities. For example, IWerner et al.l (|2008| ) find that 
most masses of baryons in the A222/A223 complex are 
locked within the halos, with WHIM contributing 10%- 
20%. 

At this point we would like to emphasize that the def- 
inition of densities adopted in WHIM studies is differ- 
ent from that one used here. Our results constrain the 
WHIM aspect within the definition of density fields typ - 
ical to spectroscopic surveys (e.g., iKovac et al.l 120101 ). 
In contrast the computation of the local density of the 
Lyman Q absorbers can only be derived from their HI col- 
umn densities, resulting in a non-t rivial role played b 
the absorber's size and shape. In iPenton et al.l (|200- 
the first comparison to LSS density has been provided, 
showing the half of the absorbers reside in voids. 

Figure H] limits the importance of the WHIM compo- 
nent in regions of high density (as defined on a 10 Mpc 
scale) which is not associated with groups to 30%. There- 
fore, the situation is much more favorable toward de- 
tecting missing baryons in underdense regions, but there 



6 



the temperature of the gas will be quite low. In fact 
the major success in resolving the missing baryons has 
been achieved usin g techniques looking for colder gas 
(jPenton et al.ll2004D . 

4. CUMULATIVE HALO MASS FUNCTIONS 

Up to this point, we have used the particle distribu- 
tion of the halos to split their mass among the cells they 
occupy. In the following we use the halo position to at- 
tribute the entire halo mass to one cell. This reduces 
the accuracy of the determination of the mass fraction 
slightly (< 5%). However, that way the treatment of ha- 
los and model galaxies which are examined hereafter and 
have unknown spatial extent can be matched. 

Plenty of studies have been devoted to investigate the 
mass function of halos in the cosmological context. Early 
analytical approaches suggested that the mass function is 
universal, i.e. its shape is independent of time and back- 
ground cosmo logy if adequate variables are used. In a 
seminal study, iPress fc' Schcchtcrl (|1974j) used the equiv- 
alent peak height, v = S c /a (see also. | Bond et all 1199 li 
iLee fc Shand"arh] H998t iSheth fc Tormenl 119990 . Such 
models can help gain insight into the statistics of grav- 
itational collapse. However, these highly nonlinear pro- 
cesses are complex, and a final validation of the mod- 
els can only be obtained by direct comparison to nu- 
merical simulations . B ased on TV-body s imulations, 
I Jenkins et all (|2001l ) and lEvrard et all (|2002t) presented 
fitting formulae for the differential mass functions accu- 
rate to ~ 10% — 20%. These studies supported the view 
that mass functi ons are indeed univ ersal. Also, the re- 
sults obtained by iLukic et al.l (|2007|) are consistent with 
a universal mass function; however, they re port a mild 
redsh if t dependence at low r edshifts (see also lReed et al.l 
120031) . IWarren et all (|2006l) found a fitting formula ac- 
curate to ~ 5% for a fixed cosmology at z = 0. Recently, 
however, based on a large set of TV -body simula t ions w ith 
different cosmological parameters iTinker et al.l (|2008|) re- 
vealed that the mass function can not be represented by a 
universal mass function at this level of accuracy. In par- 
ticular they found that the amplitude of the mass func- 
tion decreases monotonically by w 20% — 50% from z = 
to z = 2.5. Several st udies explore the mass functions 
at hi gh redshifts (e.g., iReed et all 120071 : ICohn fe White] 
l2008f ) but they do not particularly focus on universality. 

In general, analytical and numerical studies preferen- 
tially discuss the (unconditional) differential mass func- 
tion. For our purposes the integrated form, i.e. the un- 
conditional cumulative halos mass function (UHMF) is 
more useful since it directly gives the fraction of mass in 
halos relative to the total mass. Additionally, the func- 
tional forms of the unconditional differential mass func- 
tions, as given in the literature, are complex and their in- 
tegrals are even more so. Thus, for our purposes it seems 
to be a better strategy to directly examine the cumula- 
tive mass function and derive a simple fitting function 
for it. Consequently, in the first part of this section we 
introduce a fitting function for the UHMFs. In the sec- 
ond part, this function is modified to be applicable for 
the cumulative mass functions at different background 
densities, hereafter referred to as conditional cumulative 
halo mass functions (CHMFs). 

4.1. Unconditional Cumulative Halo Mass Functions 
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Fig. 6. — Upper panel: cumulative mass functions for redshifts 
z = 0.00, 0.14, 0.36, 0.69, 1.17, 1.91, and 3.06 shown as solid, dot- 
ted, dashed, dot-dashed, three-dot-dashed, and solid lines, respec- 
tively. Halo mass is given parametrically through the equivalent 
peak height v = S c /cr. At the top the corresponding mass scale for 
2 = is indicated. The red line presents a fit for the mass function 
at z = 0. The fitting function including the parameters is quoted 
in the bottom line. Lower panel: residuals of the mass functions 
and the fit at z = 0. The red line displays the residual for the mass 
function at z = which has been used for fitting. An accuracy of 
< 5% is achieved for halo, masses between 10}° and 10. 15 /i^Mq. 

Under the assumption that the initial density field is 
a Gaussian rando m field and the validity o f the spheri- 
cal collapse model, IPress fe Schechteil (|1974j) derived the 
mass fraction locked in halos above a given mass to the 
total mass in the universe, 

F PS (z/)=eric(-^, (1) 

where v is the equivalent peak height as discussed 
in $ 13.21 which, for a given cosmology, can be 
uniquely transformed into a corresponding halo mass, M . 
F{v) is equivalent to th e UHMF as introduced above. 
IPress fe Schechterl (|1974l ) argued that this parameteriza- 
tion makes the mass function universal, i.e. independent 
of evolutionary changes and cosmological parameters 
which are covered by the time and cosmology dependence 
of v. Over the last 20 years, high-resolution N-body sim- 
ulations demonstrated that Equation [T] reproduces the 
numerical halo mass function qualitatively. However, 
certain systemat ic deviations beca me apparent (e.g., 
Sheth fc Tormenl I2002t IWarren et all 120061 : [Tinker et all 
20081 ). Compared to TV-body results, Equation Q] pro- 
duces too many halos with masses corresponding to v = 1 
and to few at the high mass end. Nevertheless, we use 
the functional form of Equation Q] as template for our 
fitting formula: 

F m (is) = A eric (a v b ) , (2) 

with the fitting parameters A, a, and b. The upper panel 
of Figure [6] displays UHMFs for FoF halos derived from 
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Fig. 7. — Dependence of the halo mass fraction on environment for the redshifts z = 0, 1, and 3. Lines styles correspond the CHMFs in 
regions confined by the density contrasts, 5gio as indicated. Green lines represent a fits to the UHMF using Eq. [2] The fitting param eters 
are displayed at the bottom of the upper panels. Red lines show fits for the CHMFs using the fitting procedure described in Section 14.21 
The only difference between the upper and lower panels is the mass normalization. In the upper panels, F gives the halo mass fraction 
with respect to the total mass in the box and in the lower panels with respect to the total mass enclosed in the overdense regions. 
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Fig. 8. — Coordinate transformation used to convert uncon- 
ditional into conditional mass functions for the indicated density 
contrasts. To guide the eyes the green line for Sgw = displays 
the identity. For higher contrasts, this transformation keeps the 
mass function for large values of ln(^) unchanged yet it causes the 
mass function at low ln(i') to level off. 



the MS. Different line styles represent UHMFs derived 
from snapshots at redshift z = 0.00, 0.14, 0.36, 0.69, 1.17, 
1.91, and 3.06. The red line shows a fit to the UHMF at 
z = 0. The corresponding fitting parameters are given in 
the panel. For the computation of the ha lo masses, we 
applie d the correction formula proposed bv lWarren et al.l 
(2006?): ^corrected = N(l - N os ), where N denotes the 
number of particles within a given FoF halo. This cor- 
rection is most effective at the low mass end. The fit- 
ting range is confined by FoF halos with masses between 
2 x 10 10 H^Mq (200 particles) and 10 15 H^Mq. 
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Fig. 9. — Environmental dependence of the cumulative galaxy 
mass function. Except for the usage of galaxy masses instead of 
halo masses the plot is equal to the upper left panel in Figure [7] 
the line styles are adopted from there as well. 

The lower panel of Figure [S] shows the residuals be- 
tween the UHMFs at the given redshifts and the fit based 
on the z = mass function. The red line highlights the 
residual between the fit and the z = mass function. At 
z = and for masses in the range between 10 10 /t _1 M 
and 10 15 h~ 1 M & , the fit is accurate to the 5% level. How- 
ever, we note an increasing offset with redshift which re- 
sults in a deviation of 15% at z = 3 . 06. S imilar findings 
have been reported in iTinker et al.l (|2008f ) . The expres- 
sion, Equation H3 is an excellent fitting function for the 
cumulative mass function at z = 0. It is "universal", i.e., 
independent of cosmology at a degree to which the pa- 
rameterization by the equivalent peak height, v, permits. 





Fig. 11. — Differential galaxy mass function in different environments (black lines) separating the contribution of central (red lines) and 
non-central galaxies (green lines) in groups above the masses indicated in the upper right corner of the panels. The green lines show the 
contribution of galaxies not belonging to those groups. The particular role of central galaxies is that they provide a hump-like contribution 
to the mass function and explain a variation of that hump with the environment. 



4.2. Conditional Cumulative Halo Mass Function 

The term "conditional mass function" has been used 
to address two related problems: (1) to describe the mass 
distribution of progenitor halos which end up in a given 
halo at later times; (2) to refer to the mass distribution of 
halos within a regi o n of a given background density (cf. , 
Mo fc White! [19961 IShithl [19981 ISheth fc Lem1^ll99i 
Sheth fc Tormenll2002f ) . Here we will investigate the lat- 
ter. According to excursion set theory, the basic equation 
to approach such problems is 

F(M\S , ao) = erfc ( ± 6 -^—^.) , (3) 



where S c and a are the same as used in Equation [T] and 
do and <7n are the linear density contrast and the disper- 
sion of the background density field. It gives the fraction 
of mass in collapsed halos of mass greater than M (cor- 
responding to v = 5 c /o~) in a region that has a linear 
density contrast 5p (for a comp rehensive review includ- 
ing references see. lZentii cr 2007). In the current context, 
two difficulties arise when this equation is to be applied. 
First, the boundary surface and hence the volume of the 
overdense regions can be quite irregular which makes it 
difficult to determine ao appropriately. Second, to fa- 
cilitate the comparison with observations, we intend to 
compute the mass fraction within regions above a given 
Sn m , in other words cumulative with respect to the back- 
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ln(T(i/)) = A + ax + ( 1 -°~ a ) ln[1.0 + exp(x)] 
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TABLE 1 

CHMF FITTING PARAMETERS, A, a, b AND MASS FACTOR / m FOR 
THE BACKGROUND DENSITIES, 5g10i and REDSHIFTS 2 = 0, 1. AND 

3. 



ground density. However Equation [3] gives the halo mass 
fraction at a given So . Therefore, obtaining an expression 
suitable for our purposes would require an integration of 
Equation [3] with respect to Sq. This integration leads to 
a lengthy expression, which is too complex to be used as 
a model for a fitting formula (as done for the UHMF in 
Section B~Tj). 

Therefore, we choose a more phenomcnological ap- 
proach starting with the inspection of Figure [7l Black 
lines display UHMFs and CHMFs in regions above a 
given density contrast at redshifts of z = 0, z = 1, and 
z = 3. The CHMFs are determined by summing up the 
mass of halos which reside in cells above a given contrast 
Sqiq. The difference between the upper and lower panels 
is only in normalization. In the upper panels, we use the 
total mass in the box whereas in the lower panels we use 
the total mass within the cells above Saw- Therefore, 
the CHMFs in the upper panels lie systematically below 
the UHMFs; this is because a fraction of the total vol- 
ume is excluded by the density criterion - and so are the 
halos in it - but the mass fraction is still computed with 
respect to the total mass in the box. In the lower panels, 
the CHMFs lie systematically above the UHMFs which 
is a result of the decreasing total mass in the volume 
above the given background density. The latter is more 
physical but the former illustrates the fact that at the 
high-mass end all CHMFs display the same behavior as 
the UHMF. The green lines show fits to the UHMF. The 
values of the fitting parameters are given at the bottom 
of the upper panels (here the UHMFs are fitted sepa- 
rately for each redshift). The red lines are fits to the 
CHMFs. 

Guided by the behavior seen in the upper panels of Fig- 
ure[71 we conceive a fitting formula for the CHMFs which 
is based on a parameter-dependent coordinate transfor- 
mation, T(y), such that for an adequate set of param- 
eters F(T(u)) (here F denotes the UHMF) matches a 
given CHMF. We use the following functional form for 
these transformations, constructed to map high v values 



where x = b (m(^) — A). The values for the three pa- 
rameters A, a, and b are listed in Table [T] with the cor- 
responding fits shown as red lines in Figure [7] Figure [5] 
illustrates the coordinate transformations at z = 0. A 
determines where the deviation from one-to-one corre- 
spondence takes place, a gives the left hand side slope 
and b defines the smoothness of the transition between 
the two slopes. So far, we have described how to fit the 
CHMFs in the upper panels of Figure [7] To get the 
mass normalization right, i.e., to transform the fits of 
the upper panels into the those of the lower panels, these 
F{T(v)) have to be multiplied by the mass factor /m 
which is listed in the rightmost column of Table [TJ It 
gives the ratio of the total mass in the box to the mass 
confined to the considered overdense regions. 

Our findings may help to illustrate some of the main 
mechanisms shaping the halo mass functions in differ- 
ent environments. The top heavy shape of the CHMFs 
is caused by a relative lack of small mass halos at high 
background densities. For our way of parameterizing the 
background density, namely, by using isodensity surfaces 
to indicate all the volume interior to it, the shape of the 
cumulative mass functions at the high mass end is inde- 
pendent of environment. The only difference is induced 
by the normalization which reflects the mass confined to 
the high density regions. It causes the amplitude of the 
mass function to rise. 

As a concluding remark we note that the use of cumu- 
lative instead of differential mass functions proved valu- 
able to extract these results. First, because it directly 
gives the fraction of mass in halos above a given mass 
with respect to the total mass. In addition, it's shape is 
simple compared to the differential mass function which 
eases finding a suitable fitting strategy. Finally, it allows 
us to easily derive some of its basic properties by ana- 
lytically. For instance, it is obvious that the "wrongly 
normalized" CHMFs in the upper panels must coincide 
at the high mass end. Similarly easy to derive is that 
the fits in the lower panels must be confined by I even 
for extrapolations to smallest halo masses. The fact that 
this is obviously not the case for two of the z = 3 fits 
demonstrates the limits of the fitting procedure. How- 
ever, in general the fits behave well, i.e., the amplitudes 
remain < I even if extrapolated to small masses. 

5. CONDITIONAL GALAXY MASS FUNCTION 

The semi-analytical mod eling of galaxies i nc luded 
in the MS database (iDe Lucia fc Blaizol 120071 : 
iLemson fc Virgo Consortium! 120061 ) provides stellar 
masses for each model galaxy. This allows us to com- 
pute the conditional cumulative galaxy stellar mass 
function (CGMF) in exactly the same way as the 
CHMFs have been determined, namely by summing 
up the mass of all galaxies which reside in cells above 
a given contrast Saw- We set the lower limit for the 
galaxy masses to 10 9 /i _1 M0. Figure [5] shows the results 
of the CGMF normalized by the total mass in the box, 
which is the reason why all mass functions coincide 
at the high-mass end, equivalently to the behavior 
of the CHMFs. In general, the changes in CGMF 
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are qualitatively very similar to that of the CHMFs: 
becoming top heavy at high background densities. 

The qualitative similarities between the CHMFs and 
CGMFs suggest that the dependence of the model galaxy 
mass function on environment is caused by the corre- 
sponding dependence of the halo mass function on envi- 
ronment, rather than on direct impact of environment on 
galaxy evolution (generally referred to as nurture effects). 
To test this conjecture, Figure HU1 shows the model galaxy 
mass functions at different background densities with the 
additional restriction that they reside in dark matter ha- 
los of a given mass range as indicated by the labels right 
on top. The upper panels show the cumulative galaxy 
mass functions and the lower panels depict the GMF in 
differential form. The line styles correspond to differ- 
ent density thresholds listed in the figure. The integer 
numbers indicate the number of halos in those regions. 
Evidently, these numbers show a strong dependence on 
environment assuring that the galaxy mass functions dis- 
played in a single panel are based on very different sets of 
host halos. Nevertheless, the resulting mass functions de- 
viate by less than 10% for galaxy masses < 10 11 ft. _1 M Q . 
Therefore, the large changes in the CGMFs can almost 
entirely be accounted for by changes in the host halo 
population. Nurture effects may have a minor impact 
on shaping the model galaxy mass function in different 
environments. 

There are some shortcomings in the semianalytical 
mode l of galaxy formation used here (jDe Lucia fc Blaizotl 
120071 ). The two most important are: (1) cooling is in- 
stantaneously shut down for galaxies whose halo enters 
a larger one; (2) tidal forces are not allowed to strip off 
stars reducing the luminosity of a given galaxy. Never- 
theless, these processes do mostly affect satellite galaxies. 
Thus, we believe that the behavior seen at the high mass 
end of the CGMFs should be a reliable prediction from 
the semianalytical model. 

To provide a further illustration, we display the dif- 
ferential galaxy mass function as a function of environ- 
ment separating the contribution of central and non- 
central galaxies in groups in Figure 111! The particu- 
lar role of central galaxies is that they provide a hump- 
like contribution to the mass function and illustrate the 
variation of that hump with the environment, there- 
fore explain another im portant observational result (e.g., 
iBolzonella et al.l l2009f ). Explaining the hump- like fea- 
ture as being due to the contribution of central galaxies 
(and hence a bimodal GMF with centrals and satellites as 
its two constituent s) is somewhat d i fferen t than the ex- 
planation given by IBolzonella et al.l (|2009f ) , namely that 
the hump is due to the contributions of red galaxies to 
the total GMF. In our view, the hump is a consequence 
of halo assembly, not of the transformation of blue into 
red galaxies, although the two processes might be linked 
in the sense that the assembly of central galaxies might 
also lead to quenching of their star formation in some 
halos. We therefore predict that the hump will also be 
seen in the blue galaxy mass function alone, not only in 
the total GMF, since many halos host giant (blue) spiral 
galaxies as their central galaxy. Indeed, this bimodality 
in the blue galaxy mass function has been observed by 



iDrorv et al.l (|2009fl . 

The thermodynamical state of the accreting gas is ex- 
pected to vary with redshift, enabling a cold accretion at 
z > 2 and subsequent star formation i n the central galax y 
even in halos as massive as 1O 13 M (jDekel et al.ll2009t) . 
These are expected to be the highest peaks of density 
field at those redshifts and therefore obey our predictions 
for the behavior in high density environments. Observa- 
tional search for the transitional halo mass between the 
cold and the hot accretion mode is difficult, as this mass 
is well below the sensitivity of both spectroscopic and X- 
ray surveys for defining the galaxy groups. Instead, here 
we propose to use the shape of the galaxy mass function 
to determine at which halo mass the transition occurs. 
Since the central galaxies of groups make a hump-like 
contribution to the galaxy mass function, the blue central 
galaxies should constitute a blue hump. The location of 
the blue hump in the galaxy mass function can therefore 
be used through a comparison to numerical simulations 
to determine the transition halo mass scale for shutting 
down the star formation or environmental dependence 
of cold accretion mode in galaxy formation. In contrast 
to clustering studies, the proposed method does not in- 
duce a requirement on the data to be representative, and 
can therefore be applied to a field of any size or even a 
selected object, like a high-rcdshift supcrclustcr. 

6. CONCLUSIONS 

Using the Millennium Simulation and current schemes 
for density field reconstruction, we parameterize the con- 
ditional halo mass function. As an application, we con- 
sider the role of halos in explaining the missing baryon 
problem and the environmental dependence of galaxy 
mass functions. We show that in high-density environ- 
ments galaxy groups provide a major contribution to to- 
tal matter content. We discuss the implication of this 
result for search of missing baryons. Under the well- 
justified assumption that baryons follow dark matter, we 
show that its amount can be constrained using the obser- 
vations of galaxy groups. We also point out that the en- 
vironmental changes in galaxy mass functions are caused 
by changes in mass function of groups and, in particu- 
lar, that a hump-like features in galaxy mass function is 
produced by the central galaxies of groups. 
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